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Abstract 

We give a pedagogical introduction into the functional renormalization group treatment of disordered systems. After 
a review of its phenomenology, we show why in the context of disordered systems a functional renormalization group 
treatment is necessary, contrary to pure systems, where renormalization of a single coupling constant is sufficient. 
This leads to a disorder distribution, which after a finite renormalization becomes non-analytic, thus overcoming the 
predictions of the seemingly exact dimensional reduction. We discuss, how the non-analyticity can be measured in a 
simulation or experiment. We then construct a renormalizable field theory beyond leading order. We discuss an elastic 
manifold embedded in N dimensions, and give the exact solution for N — > oo. This is compared to predictions of 
the Gaussian replica variational ansatz, using replica symmetry breaking. We further consider random field magnets, 
and supersymmetry. We finally discuss depinning, both isotropic and anisotropic, and universal scaling function. 
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1 Introduction 

Statistical mechanics is by now a rather mature branch of physics. For pure systems like a ferromagnet, it allows 
to calculate so precise details as the behavior of the specific heat on approaching the Curie-point. We know that it 
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Figure 1: An Ising magnet at low temperatures forms a domain wall described by a function u(x) (right). An experiment on a 
thin Cobalt film (left) [6 |; with kind permission of the authors. 



diverges as a function of the distance in temperature to the Curie-temperature, we know that this divergence has the 
form of a power-law, we can calculate the exponent, and we can do this with at least 3 digits of accuracy. Best of all, 
these findings are in excellent agreement with the most precise experiments. This is a true success story of statistical 
mechanics. On the other hand, in nature no system is really pure, i.e. without at least some disorder ("dirt"). As 
experiments (and theory) seem to suggest, a little bit of disorder does not change the behavior much. Otherwise 
experiments on the specific heat of Helium would not so extraordinarily well confirm theoretical predictions. But 
what happens for strong disorder? By this we mean that disorder completely dominates over entropy. Then already 
the question: "What is the ground-state?" is no longer simple. This goes hand in hand with the appearance of so- 
called metastable states. States, which in energy are very close to the ground-state, but which in configuration-space 
may be far apart. Any relaxational dynamics will take an enormous time to find the correct ground-state, and may fail 
altogether, as can be seen in computer-simulations as well as in experiments. This means that our way of thinking, 
taught in the treatment of pure systems, has to be adapted to account for disorder. We will see that in contrast to 
pure systems, whose universal large-scale properties can be described by very few parameters, disordered systems 
demand the knowledge of the whole disorder-distribution function (in contrast to its first few moments). We show 
how universality nevertheless emerges. 

Experimental realizations of strongly disordered systems are glasses, or more specifically spin-glasses, vortex- 
glasses, electron-glasses and structural glasses (not treated here). Furthermore random-field magnets, and last not 
least elastic systems in disorder. 

What is our current understanding of disordered systems? It is here that the success story of statistical mechanics, 
with which we started, comes to an end: Despite 30 years of research, we do not know much: There are a few exact 
solutions, there are phenomenological methods (like the droplet-model), and there is the mean-field approximation, 
involving a method called replica-symmetry breaking (RSB). This method is correct for infinitely connected systems, 
e.g. the SK-model (Sherrington Kirkpatrick model), or for systems with infinitely many components. However it is 
unclear, to which extend it applies to real physical systems, in which each degree of freedom is directly coupled only 
to a finite number of other degrees of freedom. 

Another interesting system are elastic manifolds in a random medium, which has the advantage of being ap- 
proachable by other (analytic) methods, while still retaining all the rich physics of strongly disordered systems. Here, 
we review recent advances. This review is an extended version of (H|2). For lectures on the internet see (3]|4l[5). 
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2 Physical realizations, model and observables 



Before developing the theory to treat elastic systems in a disordered environment, let us give some physical realiza- 
tions. The simplest one is an Ising magnet. Imposing boundary conditions with all spins up at the upper and all spins 
down at the lower boundary (see figure 1), at low temperatures, a domain wall separates a region with spin up from 
a region with spin down. In a pure system at temperature T = 0, this domain wall is completely flat. Disorder can 
deform the domain wall, making it eventually rough again. Two types of disorder are common: random bond (which 
on a course-grained level also represents missing spins) and random field (coupling of the spins to an external random 
magnetic field). Figure 1 shows, how the domain wall is described by a displacement field u(x). Another example 
is the contact line of water (or liquid Helium), wetting a rough substrate, see figure |2] (The elasticity is long range). 
A realization with a 2-parameter displacement field u(x) is the deformation of a vortex lattice: the position of each 
vortex is deformed from x to x + u(x). A 3 -dimensional example are charge density waves. 
All these models have in common, that they are described by a displacement field 

x G M. d — > u(x) G 1^ . (1) 

For simplicity, we set N = 1 in the following. After some initial coarse-graining, the energy T~L — 7~tc\ ~\~ TYdo 
consists out of two parts: the elastic energy 

H&[u] = J d d x^(Vu(x)) 2 (2) 

and the disorder 

HdoM = J d d xV(x,u(x)) . (3) 

In order to proceed, we need to specify the correlations of disorder. Suppose that fluctuations u in the transversal 
direction scale as 

\u(x) - u{y)] 2 ~ \x - y| 2C (4) 

with a roughness-exponent £ < 1. Starting from a disorder correlator 



V{x, u)V{x', u 1 ) = f(x - x')R{u - u) (5) 




Figure 2: A contact line for the wetting of a disordered substrate by Glycerine Q. Experimental setup (left). The disorder 
consists of randomly deposited islands of Chromium, appearing as bright spots (top right). Temporal evolution of the retreat- 
ing contact-line (bottom right). Note the different scales parallel and perpendicular to the contact-line. Pictures courtesy of 
S. Moulinet, with kind permission. 
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Figure 3: Cartoon of an elastic lattice (e.g. vortex lattice) deformed by disorder. This is described by a vector u(x). 



and performing one step in the RG-procedure, one has to rescale more in the x-direction than in the ^-direction. This 
will eventually reduce f(x — x') to a ^-distribution, whereas the structure of R(u — v!) remains visible. We therefore 
choose as our starting model 

V(x,u)V(x',u') := 5 d (x - x')R(u - v!) . (6) 

There are a couple of useful observables. We already mentioned the roughness-exponent (. The second is the 
renormalized (effective) disorder. It will turn out that we actually have to keep the whole disorder distribution function 
R{u), in contrast to keeping a few moments. Other observables are higher correlation functions or the free energy. 



3 Treatment of disorder 

Having defined our model, we can now turn to the treatment of disorder. The problem is to average not the partition- 
function, but the free energy over disorder: T = — ksT In Z. This can be achieved by the beautiful replica-trick. The 
idea is to write 

lnZ= lim- fe nln2 -l N ) = lim - (Z n - I) (7) 

and to interpret Z n as the partition-function of an n times replicated system. Averaging e~ So=i ^ over disorder 
then leads to the replica-Hamiltonian 

W M = ^E / * d *\&u*{x)f - Y, J d d xR(u a (x)-u b (x)). (8) 

o=l a,b=l 

Let us stress that one could equivalently pursue a dynamic (see section [TBI or a supersymmetric formulation (section 
UTt . We therefore should not, and in fact do not encounter, problems associated with the use of the replica-trick, as 
long as we work with a perturbative expansion in R. 



4 Flory estimates 

Four types of disorder have to be distinguished, resulting in different universality classes: 

(i) Random-Bond disorder (RB): short-range correlated potential-potential correlations, i.e. short-range correlated 
R(u). 

(ii) Random-Field disorder (RF): short-range correlated force-force correlator A (it) := —R"(u). As the name 
says, this disorder is relevant for Random-field systems, where the disorder potential is the sum over all mag- 
netic fields in say the spin-up phase. 

(iii) Generic long-range correlated disorder: R{u) ~ M~ 7 . 
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(iv) Random-Periodic disorder (RP): Relevant when the disorder couples to a phase, as e.g. in charge-density waves. 
R(u) = R(u + 1), supposing that u is periodic with period 1. 

To get an idea how large the roughness Q becomes in these situations, one compares the contributions of elastic 
energy and disorder, and demands that they scale in the same way. This estimate has first been used by Flory 
for self-avoiding polymers, and is therefore called the Flory estimate. Despite the fact that Flory estimates are 
conceptually crude, they often yield a rather good approximation. For RB this gives for an iV-component field u: 

f x (Vu) 2 ~ f x vW, or L d "V ~ L d VL~ d u- N , i.e. u ~ L< with 

,rb 4 ~ rf 

.Flory 4 + _/V 



RB 4 - d 

SFIotv = A , at • (") 



For RF it is R" that is short-ranged, and we obtain 



►RF 4 - d 



^Flory 2 + N 

For LR 



(10) 



_ 4-d 

jry — 

For RP, the amplitude of u is fixed, and thus £rp = 0. 



5 Dimensional reduction 

There is a beautiful and rather mind-boggling theorem relating disordered systems to pure systems (i.e. without 
disorder), which applies to a large class of systems, e.g. random field systems and elastic manifolds in disorder. It is 
called dimensional reduction and reads as follows! 8]: 

Theorem: A d-dimensional disordered system at zero temperature is equivalent to all orders in perturbation theory 
to a pure system in d — 2 dimensions at finite temperature. Moreover the temperature is (up to a constant) nothing 
but the width of the disorder distribution. A simple example is the 3-dimensional random-field Ising model at zero 
temperature; according to the theorem it should be equivalent to the pure 1 -dimensional Ising-model at finite temper- 
ature. But it has been shown rigorously, that the former has an ordered phase, whereas we have all solved the latter 
and we know that there is no such phase at finite temperature. So what went wrong? Let us stress that there are no 
missing diagrams or any such thing, but that the problem is more fundamental: As we will see later, the proof makes 
assumptions, which are not satisfied. Nevertheless, the above theorem remains important since it has a devastating 
consequence for all perturbative calculations in the disorder: However clever a procedure we invent, as long as we 
do a perturbative expansion, expanding the disorder in its moments, all our efforts are futile: dimensional reduction 
tells us that we get a trivial and unphysical result. Before we try to understand why this is so and how to overcome it, 
let us give one more example. Dimensional reduction allows to calculate the roughness-exponent £ defined in equa- 
tion (01). We know (this can be inferred from power-counting) that the width u of a d-dimensional manifold at finite 
temperature in the absence of disorder scales as u ~ x^ 2 ~ d ^ 2 . Making the dimensional shift implied by dimensional 
reduction leads to 

[u{x) - u{Q)] 2 ~ x 4 " d = x 2C i.e. C = • (12) 



6 The Larkin-length, and the role of temperature 

To understand the failure of dimensional reduction, let us turn to an interesting argument given by Larkin [9]. He 
considers a piece of an elastic manifold of size L. If the disorder has correlation length r, and characteristic potential 
energy /, this piece will typically see a potential energy of strength 

£D = /f-Y • (13) 
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On the other hand, there is an elastic energy, which scales like 

E el = cL d ~ 2 . (14) 
These energies are balanced at the Larkin-length L = L c with 

More important than this value is the observation that in all physically interesting dimensions d < 4, and at scales 
L > L c , the membrane is pinned by disorder; whereas on small scales the elastic energy dominates. Since the disorder 
has a lot of minima which are far apart in configurational space but close in energy (metastability), the manifold can 
be in either of these minimas, and the ground-state is no longer unique. However exactly this is assumed in e.g. the 
proof of dimensional reduction; as is most easily seen in its supersymmetric formulation, see [ 10] and section [P71 

Another important question is, what the role of temperature is. In (O, we had supposed that u scales with the 
systems size, u ~ L*\ From the first term in ((8]> we conclude that 

T~L e , 6 = d-2 + 2( (16) 

Temperature is irrelevant when 9 > 0, which is the case for d > 2, and when £ > even below. The RG fixed point 
we are looking for will thus always be at zero temperature. 

From the second term in (O we conclude that disorder scales as 

R ~ L ^-4+4C . (17) 

This is another way to see that d = 4 is the upper critical dimension. 



7 The functional renormalization group (FRG) 

Let us now discuss a way out of the dilemma: Larkin's argument (section [6]) or Eq. (fTTT ) suggests that 4 is the upper 
critical dimension. So we would like to make an e = 4 — d expansion. On the other hand, dimensional reduction tells 
us that the roughness is £ = (see (fl"2~l) ). Even though this is systematically wrong below four dimensions, it tells 
us correctly that at the critical dimension d = 4, where disorder is marginally relevant, the field u is dimensionless. 
This means that having identified any relevant or marginal perturbation (as the disorder), we find immediately another 
such perturbation by adding more powers of the field. We can thus not restrict ourselves to keeping solely the first 
moments of the disorder, but have to keep the whole disorder-distribution function R{u). Thus we need afunctional 
renormalization group treatment (FRG). Functional renormalization is an old idea going back to the seventies, and 
can e.g. be found in ifTTTl . For disordered systems, it was first proposed in 1986 by D. Fisher lH"2l . Performing an 
infinitesimal renormalization, i.e. integrating over a momentum shell a la Wilson, leads to the flow diR{u), with 
(e = 4 - d) 

d e R(u) = (e - 4C) R(u) + (uR'(u) + ^R"(u) 2 - R"(u)R"(0) . (18) 

The first two terms come from the rescaling of R in Eq. (fTTT ) and u respectively. The last two terms are the result of 
the 1-loop calculations, which are derived in appendix lAl 

More important than the form of this equation is it actual solution, sketched in figure |4] After some finite renor- 
malization, the second derivative of the disorder R"(u) acquires a cusp at u = 0; the length at which this happens is 
the Larkin-length. How does this overcome dimensional reduction? To understand this, it is interesting to study the 
flow of the second and forth moment. Taking derivatives of (PTST ) w.r.t. u and setting u to 0, we obtain 

d e R"(0) = (e - 2C) R"(0) + R"'(0) 2 — ► (e- 20^(0) (19) 
d e R""(0) = eR""(0) + Si?""^) 2 + AR"'(0)R""'{0) — > eR""{0) + 3R""{0) 2 . (20) 
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Figure 4: Change of —R"(u) under renormalization and formation of the cusp. 



Since R(u) is an even function, and moreover the microscopic disorder is smooth (after some initial averaging, if 
necessary), R"'(0) and R""'(0) are 0, which we have already indicated in Eqs. (fl9l ) and (l20l) . The above equations for 
#"(0) and R""(0) are in fact closed. Equation (HU) tells us that the flow of R"(0) is trivial and that C = e/2 = 
This is exactly the result predicted by dimensional reduction. The appearance of the cusp can be inferred from 
equation (l20l ). Its solution is 

*""(% = Lacff-i)/. > C := R ""W\e=o PI) 

Thus after a finite renormalization R""(0) becomes infinite: The cusp appears. By analyzing the solution of the 
flow-equation (fT8T >. one also finds that beyond the Larkin-lengfh R"(0) is no longer given by GU) with ^"'(O) 2 = 0. 
The correct interpretation of (fl9l ). which remains valid after the cusp-formation, is (for details see below) 

d e R"(0) = (e - 2C) R"(0) + R"'(0 + ) 2 . (22) 

Renormalization of the whole function thus overcomes dimensional reduction. The appearance of the cusp also 
explains why dimensional reduction breaks down. The simplest way to see this is by redoing the proof for elastic 
manifolds in disorder, which in the absence of disorder is a simple Gaussian theory. Terms contributing to the 2-point 
function involve R"(0), TR""(0) and higher derivatives of R(u) at u = 0, which all come with higher powers of T. 
To obtain the limit of T — > 0, one sets T = 0, and only R"(0) remains. This is the dimensional reduction result. 
However we just saw that R""(0) becomes infinite. Thus R""(0)T may also contribute, and the proof fails. 



8 Measuring the cusp 

Until now the function R(u), a quantity central to the FRG, was loosely described as an effective disorder correlator, 
which evolves under coarse-graining towards a non-analytic shape. It turns out that it can be given a precise definition 
as an observable lfT3l . Hence it can directly be computed in the numerics, as we will discuss below, and in principle, 
be measured in experiments. The cusp therefore is not a theoretical artefact, but a real property of the system, 
related to singularities or shocks, which arise in the landscape of pinning forces. Moreover, these singularities are 
unavoidable for a glass with multiple metas table states. 

Consider our interface in a random potential, and add an external quadratic potential well, centered around w: 

2 

nZM = ^Hx) - wf + n el [ u ] + HdoM • (23) 

In each sample (i.e. disorder configuration), and given w, one finds the minimum energy configuration. This ground 
state energy is 

V>):=min?Ct[«]. ( 24 ) 

u(x) 

It varies with w as well as from sample to sample. Its second cumulant 

V{w)V{w>) = L d R(w - w') (25) 
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z (z/4 for RB) 



Figure 5: Filled symbols show numerical results for Y(z), a normalized form of the interface displacement correlator —R"(u) 
[Eq. (126*1)1. for D = 2 + 1 random field (RF) and D = 3 + 1 random bond (RB) disorders. These suggest a linear cusp. The 
inset plots the numerical derivative Y'{z), with intercept F'(0) w —0.807 from a quadratic fit (dashed line). Open symbols plot 
the cross-correlator ratio Y s (z) = A 12(2) / Au(0) between two related copies of RF disorder. It does not exhibit a cusp. The 
points are for confining wells with width given by M 2 = 0.02. Comparisons to 1-loop FRG predictions (curves) are made with 
no adjustable parameters. Reprinted from |fl5l . 



defines a function R(w) which is proven [13] to be the same function computed in the field theory, defined from the 
zero-momentum action [14]. 

Physically, the role of the well is to forbid the interface to wander off to infinity. The limit of small m is then 
taken to reach the universal limit. The factor of volume L d is necessary, since the width u 2 of the interface in the 
well cannot grow much more than m^^. This means that the interface is made of roughly L/L m pieces of internal 
size L m m pinned independently: (1251) hence expresses the central limit theorem and R(w) measures the second 
cumulant of the disorder seen by any one of the independent pieces. 

The nice thing about (T25T ) is that it can be measured. One varies w and computes (numerically) the new ground- 
state energy; finallying averaging over many realizations. This has been performed recently in lfT3Tl using a powerful 
exact-minimization algorithm, which finds the ground state in a time polynomial in the system size. In fact, what was 
measured there are the fluctuations of the center of mass of the interface u(w) = L~ d f d d x uo(x; w): 

[w - u(w)][w' - u(w')] C = m~' i L- d A(w - w') (26) 
which measures directly the correlator of the pinning force A (it) = —R"{u). To see why it is the total force, write 



1 
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1 1 

IF, D=3+l, L=16, M 2 =0.02 " 

Sample 1 

---- Sample 2 (5=0.1) " 


1 




1 



1 



w 

Figure 6: Discontinuous positions, "shocks", in w — u w as a function of w. Reprinted from [ 15 1. 
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the equilibrium condition for the center of mass m 2 [w — u(w)] + L~ d f d d x F(x, u) = (the elastic term vanishes if 
we use periodic b.c). The result is represented in figure[5] It is most convenient to plot the function Y = A(u) /A(0) 
and normalize the u-axis to eliminate all non-universal scales. The plot in figure [5] is free of any parameter. It has 
several remarkable features. First, it clearly shows that a linear cusp exists in any dimension. Next it is very close to 
the 1-loop prediction. Even more remarkably, as detailed in 031 . the statistics is good enough to reliably compare 
the deviations to the 2-loop predictions obtained in section [TTI 

What is the physics of the cusp in A (it)? One easily sees in a zero-dimensional model, i.e. a particle on a line, 
d = 0, that as w increases, the position of the minimum u(w) increases smoothly, except at some discrete set of 
positions w = w s , where the system switches abruptly between two distant minima. Formally, one can show that the 
landscape of the force —V'(w) evolves, as the mass is lowered, according to a Burgers equation, known to develop 
finite-time singularities called "shocks". For details on this mapping see lfl6l[P3lL For an interface these shocks also 
exist, as can be seen on figure [6l Note that when we vary the position w of the center of the well, it is not a real 
motion. It just means to find the new ground state for each w. Literally "moving" w is another very interesting 
possibility, and will be discussed in section [TBI devoted to depinning lfT7l[T8l . 

9 Rounding the cusp 

As we have seen, a cusp non-analyticity necessary arises at zero temperature, due to the switch-over between many 
metastable states. Interestingly, this cusp can be rounded by several effects: By non-zero temperature T > 0, chaos, 
or a non-zero driving velocity (in the dynamics discussed below). It is easy to include the effect of temperature in the 
FRG equation to one loop IT9l : 

d t R{u) = (e - 4C) R(u) + (uR'(u) + ^R"(u) 2 - R"(u)R"(0) + f e R"(u) . (27) 

Ti = Te~ ee is the dimensionless temperature. It finally flows to zero, since temperature is an irrelevant variable as 
discussed above. Although irrelevant, it has some profound effect. Clearly the temperature in (1271 acts as a diffusive 
term smoothening the cusp. In fact, at non-zero temperature there never is a cusp, and R(u) remains analytic. The 
convergence to the fixed point is non-uniform. For u fixed, R(u) converges to the zero-temperature fixed point, 
except near u = 0, or more precisely in a boundary layer of size u ~ Tt, which shrinks to zero in the large-scale 
limit. Non-trivial consequences are: The curvature blows up as R""(0) ~ e */T ~ L e /T. One can show that this is 
related to the existence of thermal excitations ("droplets") in the statics [20] and of "barriers" in the dynamics, which 
grow as L e (2TJ. 

Another case where rounding occurs is for "disorder chaos". Disorder chaos is the possibility of a system to have 
a completely different ground state at large scales, upon a very slight change in the microscopic disorder (for instance 
changing slightly the magnetic field in a superconductor). Not all types of disorder exhibit chaos. Its presence in spin 
glasses is still debated. Recently it was investigated for elastic manifolds, using FRG [22]. One studies a model with 
two copies, i = 1, 2, each seeing slightly different disorder energies Vi(x, u(x)) in Eq. ©. The latter are mutually 
correlated gaussian random potentials with a correlation matrix 

V i (x,u)V j (x',u') = 5 d (x - x')R lj {u - v!) . (28) 

At zero temperature, the FRG equations for Ru(u) = R22 (u) are the same as in (fT8l) . The one for the cross-correlator 
Ru{u) satisfies the same equation as (|27T ) above, with Tg is replaced by T := R'^ity ~ R"i (0)- It i s some kind of 
fictitious temperature, whose flow must be determined self-consistently from the two FRG equations. As in the case 
of a real temperature, it results in a rounding of the cusp. The physics of that is apparent from figure [6l which shows 
the set of shocks in two correlated samples. Since they are slightly and randomly displaced from each other, the cusp 
is rounded. 

Chaos is obtained when T grows with scale, and occurs on scales larger than the so-called overlap length. The 
mutual correlations CV,(x — x') = ([u l (x) — u l {x')][v?{x) — uJ(x')]) behave as Cij(x) = x 2 ^f(8x a ), where 5 
quantifies the difference between the two disorders at the microscopic level. Cij(x) decays at large distance as 

dj{x) ~x 2 <-"Q 
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10 Beyond 1 loop 



Functional renormalization has successfully been applied to a bunch of problems at 1-loop order. From a field theory, 
we however demand more. Namely that it 

• allows for systematic corrections beyond 1-loop order 

• be renormalizable 

• and thus allows to make universal predictions. 

However, this has been a puzzle since 1986, and it has even been suggested that the theory is not renormalizable due 

3 " 

to the appearance of terms of order ez [23]. Why is the next order so complicated? The reason is that it involves 
terms proportional to R"'(0). A look at figure 3 explains the puzzle. Shall we use the symmetry of R(u) to conclude 
that R"'(0) is 0? Or shall we take the left-hand or right-hand derivatives, related by 

R'"{0 + ) := lim R"'{u) = - lim R"'{u) =: -i?'"((T). (29) 

u— >0 u^O 

In the following, we will present our solution of this puzzle, obtained at 2-loop order, at large N, and in the driven 
dynamics. 



11 Results at 2-loop order 

For the flow-equation at 2-loop order, the result is Ifl4l |24l |25l |26l 1271 

d t R(u) = (e - 4C) R(u) + (uR'{u) + ^R"(u) 2 - R"{u)R"(0) 

+1 _ #'(o)) R'"( u f - l -R"'(Q+ fR"(u) . (30) 

The first line is the result at 1-loop order, already given in (fl~8T ). The second line is new. The most interesting term 
is the last one, which involves R"'(0 + ) 2 and which we therefore call anomalous. The hard task is to fix the pref- 
actor (— I). We have found five different prescriptions to calculate it: The sloop-algorithm, recursive construction, 
reparametrization invariance, renormalizability, and potentiality |[24l[T4l . For lack of space, we restrain our discussion 
to the last two ones. At 2-loop order the following diagram appears 

\ [R"(u) - R"(0)) R"'(uf - l -R"{u)R'"{Q+f (31) 



leading to the anomalous term. The integral (not written here) contains a sub-divergence, which is indicated by the 
box. Renormalizability demands that its leading divergence (which is of order 1/e 2 ) be canceled by a 1-loop counter- 
term. The latter is unique thus fixing the prefactor of the anomalous term. (The idea is to take the 1-loop correction 
5R in Eq. (|8TT > and replace one of the R" in it by 5R" itself, which the reader can check to leading to the terms given 
in (I3TI ) plus terms which only involve even derivatives.) 

Another very physical demand is that the problem remain potential, i.e. that forces still derive from a potential. 
The force-force correlation function being —R"(u), this means that the flow of R'(0) has to be strictly 0. (The 
simplest way to see this is to study a periodic potential.) From (TTTb one can check that this does not remain true if 
one changes the prefactor of the last term in (fTTT> : thus fixing it. 

Let us give some results for cases of physical interest. First of all, in the case of a periodic potential, which is 
relevant for charge-density waves, the fixed-point function can be calculated analytically as (we choose period 1, the 
following is for u G [0, 1]) 

R*(u) = - ( ±- + + 0(e*)) u\l - uf + const. (32) 




10 



c 


one loop 


two loop 


estimate 


simulation and exact 


d= 3 
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0.215 ±0.01 
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d= 2 
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0.625 


0.687 


0.67 ±0.02 


2/3 



Figure 7: Results for £ in the random bond case. 



This leads to a universal amplitude. In the case of random-field disorder (short-ranged force-force correlation func- 
tion) ( = |, equivalent to the Flory estimate (flOl) . For random-bond disorder (short-ranged potential-potential 
correlation function) we have to solve (l30l) numerically, with the result 

C = 0.20829804e + 0.006858e 2 + 0(e 3 ) . (33) 

This compares well with numerical simulations, see figure [7] It is also surprisingly close, but distinctly different, 
from the Flory estimate (©, £ = e/5. 

12 Finite JV 

Up to now, we have studied the functional RG for one component N = 1. The general case of finite iV is more 
difficult to handle, since derivatives of the renormalized disorder now depend on the direction, in which this derivative 
is taken. Define amplitude u := \u\ and direction u := u/\u\ of the field. Then deriving the latter variable leads 
to terms proportional to 1/u, which are diverging in the limit of u — > 0. This poses additional problems in the 
calculation, and it is a priori not clear that the theory at N ^ 1 exists, supposed this is the case for N = 1. At 1-loop 
order everything is well-defined [23]. We have found a consistent RG-equation at 2-loop order |29l : 



d e R(u) = (e - AQR(u) + (uR'(u) + -R"(u) 2 - R"(0)R"(u) + - 1 R ' {u) ( - 2R"(0)) 

2 2 u \ u J 



^ (&<( \ R 'Vn^) &»( ^ i N-l (R'( U )- U R"( U )) 2 (2R'( U )+u(R"( U )-3R"(0))) 
-R"'(0 + 



— - — R (u) H 



(34) 



The first line is the 1-loop equation, given in [23]. The second and third line represent the 2-loop equation, with the 
new anomalous terms proportional to R'"(0 + ) 2 (third line). 




Figure 8: Results for the roughness £ at 1- and 2-loop order, as a function of the number of components N. 
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The fixed-point equation (l34l can be integrated numerically, order by order in e. The result, specialized to directed 
polymers, i.e. e = 3 is plotted on figure [8] We see that the 2-loop corrections are rather big at large N, so some doubt 
on the applicability of the latter down to e = 3 is advised. However both 1- and 2-loop results reproduce well the 
two known points on the curve: £ = 2/3 for N = I and £ = for N = oo. The latter result will be given in section 
[T3l Via the equivalence ll30l of the directed-polymer problem in N dimensions treated here and the KPZ-equation 
of non-linear surface growth in N dimensions, which relate the roughness exponent ( of the directed polymer to the 
dynamic exponent zkpz i n tne KPZ-equation via ( = x , we know that ((N = 1) = 2/3. 

The line £ = 1/2 given on figure [8] plays a special role: In the presence of thermal fluctuations, we expect the 
roughness-exponent of the directed polymer to be bounded by ( > 1/2. In the KPZ-equation, this corresponds to 
a dynamic exponent zkpz = 2, which via the exact scaling relation zrpz + Ckpz = 2 is an upper bound in the 
strong-coupling phase. The above data thus strongly suggest that there exists an upper critical dimension in the KPZ- 
problem, with d uc 2.4. Even though the latter value might be an underestimation, it is hard to imagine what can go 
wrong qualitatively with this scenario. The strongest objections will probably arise from numerical simulations, such 
as EH- However the latter use a discrete RSOS model, and the exponents are measured for interfaces, which in large 
dimensions have the thickness of the discretization size, suggesting that the data are far from the asymptotic regime. 
We thus strongly encourage better numerical simulations on a continuous model, in order to settle this issue. 

13 Large N 

In the last sections, we have discussed renormalization in a loop expansion, i.e. expansion in e. In order to indepen- 
dently check consistency, it is good to have a non-perturbative approach. This is achieved by the large- N limit, which 
can be solved analytically and to which we turn now. We start from 

T-t[u,f\ = — ^ / u a {x) (-V 2 +m 2 ) u a (x) - ^21 Ja{x)u a (x) 

Z± a=l Jx a=l Jx 
n 

-^2 E B{{u a (x)-u b (x)f) . (35) 

Z± a,b=l Jx 

where in contrast to d8), we use an iV-component field u. For N = 1, we identify B(u 2 ) = R(u). We also have 
added a mass m to regularize the theory in the infra-red and a source j to calculate the effective action T(n) via a 
Legendre transform. For large N the saddle-point equation reads |f32l 

B\u 2 ab ) = B' (u 2 ab + 2TI% + 4h[B'(u 2 ab ) - B'(0)]) . (36) 

This equation gives the derivative of the effective (renormalized) disorder B as a function of the (constant) background 
field u 2 b = (u a — u b ) 2 in terms of: the derivative of the microscopic (bare) disorder B, the temperature T and the 
integrals I n := J k {k ^ m , f . 

The saddle-point equation can again be turned into a closed functional renormalization group equation for B by 
taking the derivative w.r.t. m: 

d,B{x) = ~B{x) = (e - AQB{x) + 2(xB'(x) + \&{x) 2 - B'(x)&(0) + (37) 
dm 2 e + B"(0) 

This is a complicated nonlinear partial differential equation. It is therefore surprising, that one can find an analytic 
solution. (The trick is to write down the flow-equation for the inverse function of B'(x), which is linear.) Let us only 
give the results of this analytic solution: First of all, for long-range correlated disorder of the form B'{x) ~ x~ 7 , 
the exponent £ can be calculated analytically as £ = 2 (i+ 7 ) • ^ a g rees with the replica-treatment in |[33l . the 1-loop 
treatment in 11231 . and the Flory estimate (fTT|) . For short-range correlated disorder, £ = 0. Second, it demonstrates 
that before the Larkin-length, B(x) is analytic and thus dimensional reduction holds. Beyond the Larkin length, 
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B"(0) = oo, a cusp appears and dimensional reduction is incorrect. This shows again that the cusp is not an artifact 
of the perturbative expansion, but an important property even of the exact solution of the problem (here in the limit 
of large N). 



14 Relation to Replica Symmetry Breaking (RSB) 



There is another treatment of the limit of large N given by Mezard and Parisi [33). They start from d35T ) but without 
a source-term j. In the limit of large N, a Gaussian variational ansatz of the form 



n s [u] = / u a {x) (-V 2 +m 2 ) u a (x) -—^2 a ab u a {x)u b (x) 

a=l Jx a,b=l 



(38) 



becomes exact. The art is to make an appropriate ansatz for a ab . The simplest possibility, a ab = a for all a ^ b 
reproduces the dimensional reduction result, which breaks down at the Larkin length. Beyond that scale, a replica 
symmetry broken (RSB) ansatz for a ab is suggestive. To this aim, one can break a ab into four blocks of equal size, 
choose one (variationally optimized) value for the both outer diagonal blocks, and then iterate the procedure on the 
diagonal blocks, resulting in 

(\ 




(39) 



One finds that the more often one iterates, the better the result becomes. In fact, one has to repeat this procedure 
infinite many times. This seems like a hopeless endeavor, but Parisi has shown that the infinitely often replica 
symmetry broken matrix can be parameterized by a function [a](z) with z G [0, 1]. In the SK-model, z has the 
interpretation of an overlap between replicas. While there is no such simple interpretation for the model (l38l) . we 
retain that z = describes distant states, whereas z = 1 describes nearby states. The solution of the large- N 
saddle-point equations leads to the curve depicted in figure 6. Knowing it, the 2-point function is given by 



(UkU-k) 



1 



k 2 + 



1 + 



1 

o z 2 k 2 + [a] (z) + m 2 



(40) 



[a](z.) + m 



from UV-cutoff 




IR-cutoff 



z, n z c l 

Figure 9: The function [a] (u) + m 2 as given in 
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The important question is: What is the relation between the two approaches, which both declare to calculate the same 
2-point function? Comparing the analytical solutions, we find that the 2-point function given by FRG is the same as 
that of RSB, if in the latter expression we only take into account the contribution from the most distant states, i.e. those 
for z between and z m (see figure [9]>. To understand why this is so, we have to remember that the two calculations 
were done under quite different assumptions: In contrast to the RSB-calculation, the FRG-approach calculated the 
partition function in presence of an external field j, which was then used to give via a Legendre transformation the 
effective action. Even if the field j is finally tuned to 0, the system will remember its preparation, as is the case 
for a magnet: Preparing the system in presence of a magnetic field will result in a magnetization which aligns with 
this field. The magnetization will remain, even if finally the field is turned off. The same phenomena happens here: 
By explicitly breaking the replica-symmetry through an applied field, all replicas will settle in distant states, and the 
close states from the Parisi-function [a] (z) + m 2 (which describes spontaneous RSB) will not contribute. However, 
we found that the full RSB -result can be reconstructed by remarking that the part of the curve between z m and z c is 
independent of the infrared cutoff m, and then integrating over m P2l (m c is the mass corresponding to z c ): 



(u k u. 



RSB 
fc=0 



R' m (0) 



+ 



mi 



1 



(41) 



We also note that a similar effective action has been proposed in |fl6l . While it agrees qualitatively, it does not 
reproduce the correct FRG 2-point function, as it should. 

15 Corrections at order 1/N 



In a graphical notation, we find (34 

SB® -- 



o 




(42) 



(43) 



where the explicit expressions are given in 13411 . 

By varying the IR-regulator, one can derive a /3-function at order 1/N, see [34]. At T = 0, it is UV-convergent, 
and should allow to find a fixed point. We have been able to do this at order e, showing consistency with the 1-loop 
result, see section [12] Other dimensions are more complicated. 

A /3-function can also be defined at finite T. However since temperature is an irrelevant variable, it makes the 
theory non-renormalizable, i.e. in order to define it, one must keep an explicit infrared cutoff. These problems have 
not yet been settled. 



16 Depinning transition 

Another important class of phenomena for elastic manifolds in disorder is the so-called "depinning transition": Ap- 
plying a constant force to the elastic manifold, e.g. a constant magnetic field to the ferromagnet mentioned in the 
introduction, the latter will only move, if a certain critical threshold force f c is surpassed, see figure [lOl (This is 
fortunate, since otherwise the magnetic domain walls in the hard-disc drive onto which this article is stored would 
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move with the effect of deleting all information, depriving the reader from his reading.) At / = f c , the so-called 
depinning transition, the manifold has a distinctly different roughness exponent £ (see Eq. ((Hi) from the equilibrium 
(/ = 0). For / > f c , the manifold moves, and close to the transition, new observables and corresponding exponents 
appear: 

• the dynamic exponent z relating correlation functions in spatial and temporal direction 

t ~ x z 

• a correlation length £ set by the distance to f c 

a~\f-f c r 

• furthermore, the new exponents are not all independent, but satisfy the following exponent relations 051 

P = v{z - C) v = (44) 

The equation describing the movement of the interface is 

d t u(x,t) = (V 2 + m 2 )u(x, t) + F(x, u(x, t)) , F(x,u) = —d u V(x,u) (45) 

This model has been treated at 1-loop order by Natterman et al. 051 and by Narayan and Fisher (36]. The 1-loop 
flow-equations are identical to those of the statics. This is surprising, since physically, the phenomena at equilibrium 
and at depinning are quite different. There is even the claim by [36], that the roughness exponent in the random 
field universality class is exactly £ = e/3, as it is in the equilibrium random field class. After a long debate among 
numerical physicists, the issue is today resolved: The roughness is significantly larger, and reads e.g. for the driven 
polymer £ = 1.25, instead of £ = 1 as predicted in ll36l . Clearly, a 2-loop analysis 071 is necessary, to resolve these 
issues. Such a treatment starts from the dynamic action 

S=l u(x,t)(d t -V 2 + m 2 )u(x,t) + u(x,t)A(u(x,t) - u(x,t'))u(x,t') , (46) 

J x,t Jx,t,t' 

where the "response field" u(x, t) enforces the equation of motion (1431 ) and 

F{x,u)F(x',u') = A(u - u')5 d (x - x') = -R"{u - u')5 d {x - x') (47) 

is the force-force correlator, leading to the second term in (l46l ). As in the statics, one encounters terms proportional 
to A'(0 + ) = — R"'(0 + ). Here the sign-problem can uniquely be solved by observing that the membrane only jumps 
ahead, 




f c f 

Figure 10: Velocity of a pinned interface as a function of the applied force. Zero force: equilibrium. / = f c : depinning. 
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Practically this means that when evaluating diagrams containing A(u(x, t) — u(x, t')), one splits them into two pieces, 
one with t < t' and one with t > t'. Both pieces are well defined, even in the limit of t — ► if. As the only tread-off of 
this method, diagrams can become complicated and difficult to evaluate; however they are always well-defined. 

Physically, this means that we approach the depinning transition from above. This is reflected in (l48l by the 
fact that u(x,t) may remain constant; and indeed correlation-functions at the depinning transition are completely 
independent of time 071 . On the other hand a theory for the approach of the depinning transition from below (/ < f c ) 
has been elusive so far. 

At the depinning transition, the 2-loop functional RG reads |[24l l37ll 



d e R(u) = (e - 4()R(u) + (uR'(u) + ~R"(u) 2 -R"(u)R"(Q) 



1 



+- [R"(u)-R"(0)]R 



»>r u f + t R >»(Q+fR»(u) 



(49) 



First of all, note that it is a priori not clear that the functional RG equation, which is a flow equation for A (it) = 
—R"{u) can be integrated to a functional RG-equation for R(u). We have chosen this representation here, in order 
to make the difference to the statics evident: The only change is in the last sign on the second line of (|49l . This has 
important consequences for the physics: First of all, the roughness exponent ( for the random-field universality class 
changes from ( = | to 



( = -(l + 0.14331e + ...) 



(50) 



Second, the random-bond universality class is unstable and always renormalizes to the random-field universality class, 
as is physically expected: Since the membrane only jumps ahead, it always experiences a new disorder configuration, 
and there is no way to know of whether this disorder can be derived from a potential or not. Generalizing the 
arguments of section [U it has recently been confirmed numerically that both RB and RF disorder flow to the RF fixed 
point lfT7l[T8l . and that this fixed point is very close to the solution of d49l , see figure[l2] This non-potentiality is most 
strikingly observed in the random periodic universality class, which is the relevant one for charge density waves. The 
fixed point for a periodic disorder of period one reads (remember A (it) = —R"(u)) 



2 / 2 ' 

A*(u) = — + — -(- + —] u(l 
y ' 36 108 \ 6 n ' v 



(51) 



Integrating over a period, we find (suppressing in F(x,u) the dependence on the coordinate x for simplicity of 
notation) 



J d«A*(u) = J\uF{u)F{u') = 



(52) 



In an equilibrium situation, this correlator would vanish, since potentiality requires Q duF(u) = 0. Here, there are 
non-trivial contributions at 2-loop order (order e 2 ), violating this condition and rendering the problem non-potential. 
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Figure 1 1 : The critical exponents at the depinning transition, for short range elasticity (left) and for long range elasticity (right). 
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Figure 12: Universal scaling form Y (z) for A(u) for RB and RF disorder. Reprinted from 1 18 1. 



This same mechanism is also responsible for the violation of the conjecture C = f > which could be proven on the 
assumption that the problem remains potential under renormalization. Let us stress that the breaking of potentiality 
under renormalization is a quite novel observation here. 

The other critical exponents mentioned above can also be calculated. The dynamical exponent z (for RF-disorder) 
reads WfiSH 

z = 2--e- 0.04321e 2 + . . . (53) 

All other exponents are related via the relation (|44l) . That the method works well even quantitatively can be inferred 
from figure [TT] 



17 Supersymmetry 

The use of n replicas in the limit n — ► to describe disordered systems is often criticized for a lack of rigor. It is 
argued that instead one should use a supersymmetric formulation. Such a formulation is indeed possible, both for the 
statics as discussed in [38], as for the dynamics, which we will discuss below. Following iflOl . one groups the field 
u(x), a bosonic auxiliary field u(x) and two Grassmanian fields ip(x) and ^(x) into a superfield U (x, 0, G): 

U(x, 9, 6) = u(x) + Qip(x) + $(z)6 + QGu{x) . (54) 

The action of the supersymmetric theory is 

Ssusy = / dQd © / U(x, 6, @){A s )U(x, 6, 6) , A s := V 2 - A(0)^^ (55) 
It is invariant under the action of the supergenerators 

What do the fields mean? Upon integrating over and G before averaging over disorder, one would obtain two 
terms, ~ f x u(x)-^^rj, i.e. the bosonic auxiliary field u(x) enforces = 0, and a second term, bilenear in ip and 

if), ~ j x "ip^^-ip which ensures that the partition function is one. (|55T > is nothing but the dimensional reduction result 
(fl2l) in super-symmetric disguise. What went wrong? Missing is the renormalization of R(u) itself, which in the 
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FRG approach leads to a flow of A(0) = — R"(0). In order to capture this, one has to look at the supersymmetric 
action of at least two copies: 

si u a] = Y / [ (i/.(i l 8 1 0)(A ( )[f a (i l e ) 8)-Jr/ f / R(u a (x,e,e)-u b ( x ,@',e')) . (5i) 

Formally, we have again introduced n replicas, but we do not take the limit of n — > 0; so criticism of the latter limit 
can not be applied here. After some slightly cumbersome calculations one reproduces the functional RG /3-function 
at 1-loop (fT8T >. (Higher orders are also possible and the SUSY method is actually helpful [38].) At the Larkin-length, 
where the functional RG produces a cusp, the flow of A(0) becomes non-trivial, given in (l22l ). Then the parameter 
A(0) in the supersymmetry generators (l56l ) is no longer a constant, and supersymmetry breaks down. This is, as was 
discussed in section (fl4l) . also the onset of replica-symmetry breaking in the gaussian variational ansatz, valid at large 
N. 

Another way to introduce a supersymmetric formulation proceeds via the super-symmetric representation of a 
stochastic equation of motion 11391 ; a method e.g. used in [40] to study spin-glasses. The action then changes to 

S[U]= I I U(x,@,Q,t)(A d )U(x,e,@,t)-l I I [ R(U(x, 0, 0, t) — U(x, &', t')) . (58) 
Jxt Je,e z Jxtt' Je,e Je',e' 

A^V> + DD, D-±, = i-e| (59) 

and is invariant under the action of the super-generators Q := and Q := -J^ + Jj, since {Q, D} = {Q, £>} = 
{Q,D} = {Q,D} = 0. Different replicas now become different times, but second cumulant still means bilocal 
in 0. However the procedure is not much different from a pure Langevin equation of motion, as in (|45T >: in the 
latter equation Ito discretization is already sufficient to ensure that the partition function is 1. The main advantage 
is the possibility to change the discretization procedure from Ito over mid-point to Stratonovich without having to 
add additional terms. In this case, supersymmetry breaking means that the system falls out of equlibrium, i.e. the 
fluctuation-dissipation theorem (which is a consequence of one of the supersymmetry generators [39]) breaks down 

ED. 



18 Random Field Magnets 

Another domain of application of the Functional RG is spin models in a random field. The model usually studied is: 

H = J d d x ^(V5) 2 + h(x) ■ S(x) , (60) 

where S(x) is a unit vector with iV-components, and S(x) 2 = 1. This is the so-called O(N) sigma model, to which 
has been added a random field, which can be taken gaussian hi(x)hj(x') = a5ijd d (x — x'). In the absence of disorder 
the model has a ferromagnetic phase for T < Tf and a paramagnetic phase above Tf. The lower critical dimension is 
d = 2 for any N > 2, meaning that below d = 2 no ordered phase exists. In d = 2 solely a paramagnetic phase exists 

for N > 2; for N = 2, the XY model, quasi long-range order exists at low temperature, with S(x)S(x') decaying as 
a power law of x — x'. 

Here we study the model directly at T = 0. The dimensional reduction theorem in section [51 which also holds 
for random field magnets, would indicate that the effect of a quenched random field in dimension d is similar to the 
one of a temperature T ~ a for a pure model in dimension d — 2. Hence one would expect a transition from a 
ferromagnetic to a disordered phase at a c as the disorder increases in any dimension d > 4, and no order at all for 
d < 4 and N > 2. This however is again incorrect, as can be seen using FRG. 

It was noticed by Fisher fiTTl that an infinity of relevant operators are generated. These operators, which corre- 
spond to an infinite set of random anisotropics, are irrelevant by naive power counting near d = 6 Il42ll43l . d = 6 is 
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the naive upper critical dimension (corresponding to d = 4 for the pure O(N) model) as indicated by dimensional 
reduction; so many earlier studies concentrated on d around 6. Because of these operators the theory is however hard 
to control there. It has been shown fl4i~l l42l l43l that it can be controlled using 1-loop FRG near d = 4 instead, which 
is the naive lower critical dimension. Recently this was extended to two loops ll44l . 

The 1-loop FRG studies directly the model with all the operators which are marginal in d = 4, of action most 
easily expressed directly in replicated form: 



d d x 



(61) 



ab 



The function R(z) parameterizes the disorder. Since the vectors are of unit norm, z = cos <p lies in the interval 
[—1,1]. One can also use the parametrization in terms of the variable <fi which is the angle between the two replicas, 
and define R(4>) = R(z = cos<ft). The original model (l60l > corresponds to R(z) ~ az. It does not remain of this 
form under RG, in fact again a cusp will develop near z = 1. The FRG flow equation has been calculated up to two 
loops, i.e. R 2 (one loop) BUSSES and R s (two loops) B^ : 



d e R(cf>) = eR{<j>) + -R"^) 2 - R"{Q)R"{<j ) ) + (N-2) 



1 R'(cf>f 



+~(R"(cf>) - R"(<S))R>"(<t>f + (N-2) 



cot i 



sin 



sm 



-R'(<j>Y 



-cot (j)R'((j))R"(0) 
5 + cos 2(p 



4 sin 



R'WR"^) 



+ 



1 



-R"(<pf 



l 



2 sin 2 <fi 4 sin (f> 

+(5 + cos 2(f)) sin 2 (j)R"(<i)) 2 



R"(0)[2(2 + cos2(j))R f ((j)) 2 - 6 sin 2<j> R' (4>)R" (4>) 



N+2 



R't'iO+fR"^) 



N-2 



cot (t>R"'(<d + ) 2 R'{$) 



2(N-2) R"(0) - R"(0) 2 + la R"'(0 + ) 2 R(<f>) 



(62) 



The constant 7 a is discussed in [44] ; the last factor proportional to R(<f>) takes into account the renormalization of 
temperature, a specific feature absent in the manifold problem. The full analysis of this equation is quite involved. 
The 1-loop part already shows interesting features. For N = 2, the fixed point was studied in |47l , and corresponds 
to the Bragg-glass phase of the XY model with quasi-long range order obtained in a d = 4 — e expansion below 
d = 4. Hence for N = 2 the lower critical dimension is d\ c < 4 and conjectured to be d\ c < 3 in [47 ]. On the other 
hand Feldman Il42ll43ll found that for N = 3, 4, . . . there is a fixed point for d = 4 + e for d > 4. This fixed point has 
exactly one unstable direction, hence was conjectured to correspond to the ferromagnetic-to-disorder transition. The 
situation at one loop is thus rather strange: For N = 2, only a stable FP which describes a unique phase exists, while 
for iV = 3 only an unstable FP exists, describing the transition between two phases. The question is: Where does the 
disordered phase go as N decreases? These results cannot be reconciled within one loop and require the full 2-loop 
analysis of the above FRG equation. 

The complete analysis [44 J shows that there is a critical value of N, N c = 2.8347408, below which the lower 
critical dimension d\ c of the quasi-ordered phase plunges below d = 4. Hence there are now two fixed points below 
d = 4. For N > N c a. ferromagnetic phase exists with lower critical dimension d\ c = 4. For N < N c one finds an 
expansion: 



d 



RF 
lc 



4 - e c « 4 - 0.1268(iV - N c ) 2 + 0((N - N c y) 



(63) 



'These results were confirmed in |45 1 (for the normal terms not proportional to R"'(0 + )) and [46 1 (with one proposition for the anomalous 
terms). 
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N>N C M = N C N<N C 




4 d 4 d d\ c 4 d 

Figure 13: Phase diagram of the RF non-linear sigma model. D = disordered, F = ferromagnetic, QLRO = quasi long-range 
order. Reprinted from [44]. 



One can also compute the exponents of the correlation function 

~ q~ 4+f > , (64) 

and once the fixed point is known, fj is given by fj = e — (N — 1)R"(0) + 3N ~ 2 R"'(0 + ) 2 . There is a similar 
exponent for the connected thermal correlation rj. Another fixed point describing magnets with random anisotropics 
(i.e. disorder coupling linearly to Si(x)Sj(x)) is studied in Il42ll43ll44ll48l . 

In this context, the existence of a quasi-ordered phase for the random-field XY model in d = 3 (the scalar version 
of the Bragg glas) has been questioned P31 . Corrections in (|63T > seem to be small and at first sight exclude the 
quasi-ordered phase in d = 3. This should however be taken with a (large) grain of salt [49]. Indeed the above 
model does not even contain topological defects (i.e. vortices) as it was directly derived in the continuum. In the 
absence of topological defects it is believed that the lower critical dimension is d\ c = 2 (with logarithmic corrections 
there). Hence the above series should converge to that value for N = 2, indicating higher order corrections to (l63l . 
Another analysis [5011 based on a FRG on the soft-spin model, which may be able to capture vortices, indicates 
d\ c (N = 2) > 3. Unfortunately, it uses a truncation of FRG which cannot be controlled perturbatively, and as a 
result, does not match the 2-loop result. It would be interesting to construct a better approximation which predicts 
accurately at which dimension the soft and hard spin model differ in their lower critical dimensions, probably when 
vortices become unbound due to disorder. 

19 More universal distributions 

As we have already seen, exponents are not the only interesting quantities: In experiments and simulations, often 
whole distributions can be measured, as e.g. the universal width distribution of an interface that we have computed at 
depinningc ll5Tll52ll . Be (u) the average position of an interface for a given disorder configuration, then the spatially 
averaged width 



X 



is a random variable, and we can try to calculate and measure its distribution P(w 2 ). The rescaled function 
defined by 

P{w 2 ) = l/w 2 ® (w 2 /vA (66) 



will be universal, i.e. independent of microscopic details and the size of the system. 



20 




0.5 1 1.5 2 2.5 3 3.5 4 



z 

Figure 14: Scaling function $>(z) for the (1 + l)-dimensional harmonic model, compared to the Gaussian approximation for 
C = 1.25. Data from gTJ. 

Supposing all correlations to be Gaussian, $>(z) can be calculated analytically. It depends on two parameters, 
the roughness exponent ( and the dimension d. Numerical simulations displayed on figure [14] show spectacular 
agreement between analytical and numerical results. As expected, the Gaussian approximation is not exact, but to 
see deviations in a simulation, about 10 5 samples have to be used. Analytically, corrections can be calculated: They 
are of order R"'(0 + ) A and small. Physically, the distribution is narrower than a Gaussian. 

There are more observables of which distributions have been calculated within FRG, or measured in simulations. 
Let us mention fluctuations of the elastic energy [53], and of the depinning force Il54ll531 . 

20 Anisotropic depinning, directed percolation, branching and all that 

We have discussed in section [TBI isotopic depinning, which as the name suggests is a situation, where the system is 
invariant under a tilt. This isotropy can be broken through an additional anharmonic elasticity 

^elastic = J ~ [V«(x)] 2 + C 4 [Vu(x)] 4 , (67) 

leading to a drastically different universality class, the so-called anisotropic depinning universality class, as found 
recently in numerical simulations ll56ll . It has been observed in simulations ll57l 1581 . that the drift- velocity of an 
interface is increased, which can be interpreted as a tilt-dependent term, leading to the equation of motion in the form 

d t u(x, t) = V 2 u(x, t) + A [Vu(x, t)} 2 + F(x, u(x, t)) + f . (68) 

However it was for a long time unclear, how this new term (proportional to A), which usually is referred to as a 
KPZ-term, is generated, especially in the limit of vanishing drift- velocity. In [59 ] we have shown that this is possible 
in a non-analytic theory, due to the diagram given in figure [T5l 




Figure 15: The diagram generating the irreversible nonlinear KPZ term with one disorder vertex and one C4 vertex (the bars 
denote spatial derivatives). 
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For anisotropic depinning, numerical simulations based on cellular automaton models which are believed to be in 
the same universality class fl60ll6Tl . indicate a roughness exponent £ 0.63 in d = 1 and £ 0.48 in d = 2. On a 
phenomenological level it has been argued ll60l loTl [62 ] that configurations at depinning can be mapped onto directed 
percolation in d = 1 + 1 dimensions, which yields indeed a roughness exponent £dp = v±/v\\ = 0.630 ± 0.001, and 
it would be intriguing to understand this from a systematic field theory. 

This theory was developed in 11591 . and we review the main results here. A strong simplification is obtained by 
going to the Cole-Hopf transformed fields 

Z(x,t):=e Xu ^ & u(x,t) = HZ [ X ' t]) . (69) 

A 

The equation of motion becomes after multiplying with AZ(x, t) (dropping the term proportional to /) 

r2 ln(Z(x,t)) 



d t Z(x, t) = V 2 Z{x, t) + XF[ x, x ) Z(x, t) (70) 



and the dynamical action (after averaging over disorder) 

v2 



S= f Z(x,t){d t -V 2 )Z(x,t)-^- f i Z(x, t)Z(x, t) A ( ln t] ln tl \ Z (x, t')Z(x, t') (71) 

This leads to the FRG flow equation at 1-loop order 

d e A{u) = (e - 2C)A(n) + £uA'(u) - A"(n) (A(«) - A(0)) - A'(n) 2 

+2AA(u)A'(0+) + 2A 2 (A(u) 2 + A(«)A(0)) (72) 

The first line is indeed equivalent to (fT8T ) using A(n) = —R"{u). The second line is new and contains the terms 
induced by the KPZ term, i.e. the term proportional to A in (I68T ). 

Equation (l72l possesses the following remarkable property: A three parameter subspace of exponential functions 
forms an exactly invariant subspace. Even more strikingly, this is true to all orders in perturbation theory [59]! The 
subspace in question is (0 < u < 1/A) 



A(u) = ^\ y a + b e~ Au + c e Au J (73) 
The FRG-flow (|72"1) closes in this subspace, leading to the simpler 3-dimensional flow: 

O e a = a + 4a 2 + 4ac + Abe (74) 
dgb = 6(1 + 6a + b + 5c) (75) 
d t c = c(l + 6a + b + 5c) (76) 

This flow has a couple of fixed points, given on figure [161 They describe different physical situations. The only 
globally attractive fixed point is SAP, describing self-avoiding polymers. This fixed point is not attainable from the 
physically relevant initial conditions, which lie (as fixed point RP) between the two separatrices given on figure [TBI 
All other fixed points are cross-over fixed points. 

In the Cole-Hopf representation, it is easy to see why the exponential manifold is preserved to all orders. Let us 
insert d73l in (TTTb . The complicated functional disorder takes a very simple polynomial form [59]. 

S= Z(x,t) (d t -V 2 ) Z(x,t) - / Z(x,t)Z{x,t')(aZ(x,t)Z(x,t')+bZ(x,t) 2 + cZ(x,t'f) . (77) 

Jxt Jx Jt<t' 

The vertices are plotted on figure [FT] It is intriguing to interpret them as particle interaction (a) and as different 
branching processes (b and c): Z destroys a particle and Z creates one. Vertex b can e.g. be interpreted as two particles 
coming together, annihilating one, except that the annihilated particle is created again in the future. However, if the 
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Figure 16: Fixed point structure for A = 2, which is a typical value. The ratio c/b is not renormalized, see d75l)-(l76ll, such that 
c/b is a parameter, fixed by the boundary conditions, especially A. The fixed points are Gaussian G, Random Periodic RP (the 
generalization of the RP fixed point for A = 0), Self-Avoiding Polymers SAP, and Unphysical U. 



annihilation process is strong enough, the reappearance of particles may not play a role, such that the interpretation 
as particle annihilation or equivalently directed percolation is indeed justified. 

One caveat is in order, namely that the fixed points described above, are all transition fixed point, and nothing 
can a priori be said about the strong coupling regime. However this is the regime seen in numerical simulations, for 
which the conjecture about the equivalence to directed percolation has been proposed. Thus albeit intriguing, the 
above theory is only the starting point for a more complete understanding of anisotropic depinning. Probably, one 
needs another level of FRG, so as standard FRG is able to treat directed polymers, or equivalently the KPZ-equation 
in the absence of disorder. 



21 Problems not treated in these notes. . . and perspectives 

Problems not treated in these notes are too numerous to list. Let us just mention some: [63] consider depinning at 
the upper critical dimension. In ifTTl . the crossover from short-ranged to long-ranged correlated disorder is treated. 
Our techniques can be applied to the statics at 3-loop order [64]. But many questions remain open. Some have 
already been raised in these notes, another is whether FRG can also be applied to other systems, as e.g. spinglasses or 
window glasses. Can FRG be used as a tool to go beyond mean-field or mode-coupling theories? Another open issue 
is the applicability of FRG beyond the elastic limit, i.e. to systems with overhangs and topological defects, non-liear 
elasticity [65], or to more general fractal curves than (directed) interfaces. For random periodic disorder in d = 2, 
temperature is marginal, and a freezing transition can be discussed (see e.g. 11661, |67l ). It would be interesting to 
connect this to methods of conformal field theory and stochastic Lowner evolution. We have to leave these problems 
for future research and as a challenge for the reader to plunge deeper into the mysteries of functional renormalization. 




Figure 17: The three vertices proportional to a, b and c in equation d77l >. 
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A Derivation of the functional RG equations 



In section [6l we had seen that 4 is the upper critical dimension. As for standard critical phenomena l39l . we now 
construct an e = (4 — d) -expansion. Taking the dimensional reduction result (PT21 in d = 4 dimensions tells us that 
the field u is dimensionless. Thus, the width a = —R"(0) of the disorder is not the only relevant coupling at small 
e, but any function of u has the same scaling dimension in the limit of e = 0, and might thus equivalently contribute. 
The natural conclusion is to follow the full function R(u) under renormalization, instead of just its second moment 
R"(0). Such an RG-treatment is most easily implemented in the replica approach: The n times replicated partition 
function becomes after averaging over disorder 



exp 



Y n 1 n \ ( 1 n 1 n f 

- ^eiM £d °w ) = ex p [~f Yl E ^ + ^Y1 1 d ° 

a=l a=l / \ a=l a, 6=1 



xR( uJx) 



u b [x. 



(78) 

Perturbation theory is constructed along the following lines (see l23l [i"4l for more details.) The bare correlation 
function, graphically depicted as a solid line, is with momentum k flowing through and replicas a and b 

T5 a b 



k 2 



The disorder vertex is 



::=/£*(«.« 

a.o 



(79) 



(80) 



The rules of the game are to find all contributions which correct R, and which survive in the limit of T = 0. At 
leading order, i.e. order R 2 , counting of factors of T shows that only the terms with one or two correlators contribute. 
On the other hand, YJ a b R(u a — Ub) has two independent sums over replicas. Thus at order R 2 four independent 
sums over replicas appear, and in order to reduce them to two, one needs at least two correlators (each contributing a 
Sab). Thus, at leading order, only diagrams with two propagators survive. These are the following (noting C{x — y) 
the Fourier transform of 1/ k 2 ): 



X 

a 



y 



-* b 




R"(u a {x) - u b {x))R"(u a {y) - u a (y))C{x - yf 



R"(u a {x) - u a (x))R"(u a (y) - u a (y))C(x - yf 



(81) 



(82) 



In a renormalization program, we are looking for the divergences of these diagrams. These divergences are localized 
at x = y, which allows to approximate R"(u a (y) — Ub(y)) by R"(u a {x) — Ub(x)). The integral f x C(x — y) 2 = 

J k jk^TnFp = ( usm § tne most convenient normalization for J k ) is the standard 1-loop diagram from </> 4 -theory. 
We have chosen to regulate it in the infrared by a mass, i.e. physically by the harmonic well introduced in section [8j 
Note that the following diagram also contains two correlators (correct counting in powers of temperature), but is 
not a 2-replica but a 3-replica sum: 

x y 

a <ZIZ> l a . (83) 

b * * c 

Taking into account combinatorial factors, the rescaling (fTTT ) of R, as well as of the field u (its dimension being 
the roughness exponent Q, we arrive at 



m^-R{u) = (e - 4QR(u) + (uR'(u) + -R"{uf - R"(u)R"{0) . 
dm 2 



(84) 



Note that the field does not get renormalized due to the exact statistical tilt symmetry u(x) — ► u(x) + ax: The bare 
action ([8]>, including the mass term, changes according to 7i 



bare 



6H\u] 



d d x 



Vu(x)a + —a 



+ m 



H h&re [u] +SH[u], with 



u(x)ax + 2 a2x2 



(85) 
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Figure 18: The toy model ([86j before (left) and after (right) the Larkin-scale. 



To render the presentation clearer, the elastic constant c set to c = 1 in equation (HI) has been introduced. The 
important observation is that all fields u involved are large-scale variables, which are also present in the renormalized 
action, changing according to TL Icn [u] — > H Icn {u] + 5H[u\. The latter can be used to define the renormalized elastic 
coefficient c ren and mass m ren . Since 5H[u] gives the change in energy both for the bare and the renormalized action 
with unchanged coefficients, c rcn = c and m ren = m, so neither elasticity nor mass changes under renormalization. 

B Why is a cusp necessary in 4 — e dimensions? 

Let us present here another simple argument due to Leon Balents ll68l . why a cusp is a physical necessity in order to 
have a e = 4 — d expansion. To this aim, consider a toy model with only one Fourier-mode u = u q 

H[u] = ^q 2 u 2 + y/eV(u) . (86) 

Since equation (fl~8l) has a fixed point of order R(u) ~ e for all e > 0, V(u) scales like yfe for e small and we have 
made this dependence explicit in d86l ) by using V(u) = y/eV{u). The only further input comes from the physics: For 
L < L c , i.e. before we reach the Larkin length, there is only one minimum, as depicted in figure [181 On the other 
hand, for L > L c , there are several minima. Thus there is at least one point for which 

d 2 

H[u] =q 2 + yfeV"{u) < . (87) 

dn z 

In the limit of e — > 0, this is possible if and only if -R""(0), which a priori should be finite for e 
infinite: 

~R""(0) = V"(u)V"(u')\ u=u , = oo . 
This argument shows that a cusp is indeed a physical necessity. 
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